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We explore the theoretical feasibility of extracting V u b from two ratios built 
from B meson inclusive partial decays, R\ = T(b — » ucs')/3T(b — > civ) and 
R2 = [T(b — ► cX) — r(6 — > cX)]/T(b — > cud'). We discuss contributions 
<D to these quantities from perturbative and nonperturbative physics, and show 
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that they can be computed with overall uncertainties at the level of 10% 



0. 



I. INTRODUCTION 



The accurate measurement of V u b is one of the most challenging theoretical and exper- 
imental problems in B physics. Its value is crucial for constraining the Unitarity Triangle 
and probing the question of whether the CKM framework is adequate for describing flavor 
physics in the standard model. The present best experimental values for V u b, from the inclu- 
sive decay B — > X U W and the exclusive process B — > pli>, are limited by model-dependence 
and other theoretical errors. New approaches to extracting V u b from inclusive and exclusive 
semileptonic decays have been proposed and are promising, but have not yet proven to be 
viable experimentally. 

In light of this situation, new methods for probing V u b are still needed. In a recent 
paper |TJ , we suggested that it would be useful to attempt to measure the inclusive production 
of "wrong sign" charm in B decays, that is, to look for evidence for the quark level transition 
b — ► ucs' . In particular, we proposed to study the ratio Ri = T(b — > ucs')/3T(b — ► civ), 
noting that the theoretical expression for this quantity is in a number of respects particularly 
well under control. (Here s' and d' are the flavor eigenstates, and we take m s = = 0.) We 
went on to compute the leading perturbative and nonperturbative corrections to the parton 
model result for R\. The analysis of Ref. [|IJ also relied implicitly on the use of parton-hadron 
duality. This assumption is common to all extractions of V u b from inclusive B decays, and 
while it is not unreasonable to expect it to hold in this case, there is no rigorous proof that 
it actually does. Perhaps the near equality of the charged and neutral B meson lifetimes 
provides some empirical evidence that duality is well respected in B decays. 

In this paper we will refine the analysis of Ref. in a number of respects. First, we 
will include complete radiative corrections to R\ at next-to-leading order, that is, all terms 
proportional to a s {mb) and ct™ +1 \ia n (mw / m>b) ■ Second, we will include a set of "enhanced" 
two loop terms, often referred to as "BLM" corrections @, which are proportional to a^flo, 
where (3$ = 11 — 2n//3 is the first coefficient in the QCD beta function. It was pointed out 
in Ref. |1J that these terms are not likely to be as large in Ri as in, for example T(B — > 
X C W), because of the cancellation of the leading renormalon ambiguity. Indeed, our explicit 
calculation will show that these terms contribute only at the level of ten percent. Third, in 
Ref. |]], we also included the leading nonperturbative contributions to the inclusive decay, 
which come from annihilation processes and are proportional to 167r 2 /g/m^. These terms 
are formally of order 1 / m| but are enhanced by the two-body, rather than three-body, phase 
space of the final state. We found that in charged B decays, these processes can contribute at 
the order of 5%, while in neutral B decays they turn out to be negligible. We will have little 
new to say about these corrections, except that we will attempt to combine the uncertainties 
from these contributions with those from the radiative corrections to obtain an overall picture 
of the reliability of the theoretical calculation. We will also take the opportunity to include 
an additional small "hybrid" contribution of order a s \n(mw/^ib)^2/fnb- 

We will also propose that it is useful to consider a second ratio, R 2 = [T(6 — > cX) — Y{b — > 
cX)\/Y(b — > cud'). We will see that R 2 is theoretically clean in a way which is similar to 
Ri, and we will extend all aspects of our analysis of -Ri to include R 2 . The experimental 
measurement of R 2 would certainly be challenging, but the challenges would be distinct 
from those that confront the measurement of R\ and this second ratio deserves separate 
consideration. 



1 



Finally, we will close by presenting an overall picture of the theoretical understanding of 
Ri and R2, with our best estimate of the remaining uncertainties and the future prospects 
for reducing them. We hope that this will provide an intriguing goal for our experimental 
colleagues to aim for. 

II. V ub FROM RATIOS OF PARTIAL DECAY RATES 

In our previous paper we proposed that the quark level process b — > ucs' would be 
a promising mode from which to extract V u b. Final states with this combination of quark 
flavors arise only from processes proportional to V u b, with no contributions from penguin 
diagrams or long distance rescattering. In the form of the ratio 

R 1 = T(b -> ucs')/3T(b -> cW) , (2.1) 

the theoretical expression is very well behaved. The phase space dependence on x c = m 2 c /m1 
is identical in the numerator and denominator, as are the leading nonperturbative corrections 
of order 1/m 2 . At tree level and in the limit rrib — > 00, then, we have simply 

Ri = \V ub /V cb \ 2 {l + 0(a s , l/m 2 b )} . (2.2) 

In Ref. UJ, we first included radiative corrections at leading logarithmic order, which has 
the effect of multiplying the expression for Ri by a factor x — 1-09. We also computed the 
leading radiative correction to the decay processes b — > ucs' and b — > civ, which although 
formally subleading is numerically substantial. The result was an expression of the form 




With x c = 0.09 and a s = a s {m h ) = 0.22, the one loop radiative correction is ga s /n = 0.21, 
so indeed it is large and should be included. However, the scale /1 at which ot s ought to be 
evaluated was not fixed by our calculation, leading to a significant remaining uncertainty. 
This can be resolved only with a full next-to-leading order calculation, which we will perform 
in the next section. We will find that the partial calculation of Ref. |]T[ was correct to within 
approximately 20%. 

Unfortunately, the experimental measurement of R\ is extremely challenging. The largest 
background to observing the quark level process b — *> ucs' is b — > ccs' , the rate for which is 
approximately a factor of 100 larger. The measurement is made more difficult by the fact 
that the only quantity which is well predicted theoretically is the ratio of fully inclusive 
rates, while many of the experimental techniques for rejecting b — > ccs' involve tagging on 
a particular hadronic final state. Although the measurement of R\ may be feasible, it will 
hardly be straightforward. Nevertheless, relevant experimental techniques already are being 
developed and the excellent capability of the BaBar and BELLE detectors to vertex 
individually the boosted B mesons also will improve the prospects for this measurement [Q . 

There is another ratio which one might consider, which avoids the necessity of rejecting 
the b — > ccs' background. Let T(b — > cX) = T(b — > cud') + T(b — ► ccs') be the fully inclusive 
production of c in nonleptonic b decays, and T(b ^ cX) = T(b — > ucs') + T(b — > ccs') be 
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the inclusive production of c. Also, define F(b — > cX') = T(b — > cud'), that is, the inclusive 
production of c without an accompanying c. Note that F(b — > ccs') < T(b — > cud'), so the 
measurement of T(6 — > cX') does not require rejecting an overwhelming background. Then 
let 

r(6^ c x') 2 v ; 

We see that T(b — > ccs') cancels in the numerator of i?2- I n terms of quark level transitions, 

r(6 -> mcsQ 

52 = r(6 -» cud') ■ (2 - 5) 

At tree level, 62 = \V u b/V c b\ 2 . Unlike in R\, there is no leading logarithmic correction to 
62, since these contribute identically to the numerator and the denominator. The leading 
radiative correction arises at order q; s (to&), and the leading nonperturbative corrections at 
order a s ln(m H //m b )A 2 / , n2^ and 16n 2 f^/ml. 

The experimental advantage of R2 is that the large b — > ccs' background cancels. The 
difficulty is that in order for R 2 to be sensitive to | V u b/V c b\ 2 , both T(b — > cX) and T(b cX') 
must be measured with an accuracy of better than 1%. This may prove to be as challenging 
as rejecting b — > ccs', or even more so, but it involves a distinct set of problems. It will 
be up to the experimental community to determine whether this measurement is feasible or 
not. 

A quantity which would be more attractive experimentally is 

r( fe _> cX) - T(b -> cX) 

R3 = r(6 - ciu) ' (2 ' 6) 

since doing so avoids the requirement of determining Y(b — > cX') precisely. Unfortunately, 
R 3 cannot be computed with the very small uncertainties of Ri and i?2- This is simply 
because, unlike R\ and R 2 , the calculation of R 3 requires that the ratio T(b cX')/Y{b — > 
civ) = Ri/82 be known with a theoretical accuracy of better than 1%. This is well beyond 
the level of precision presently attainable, due to the perturbative and nonperturbative con- 
tributions which we will discuss below. Thus, we will not consider R3 further. However, we 
do note that in anticipation of theoretical advances, it may well be useful for experimental- 
ists to measure R3 as well as Ri and R2. Alternatively, if T(b — > cX')/Y{b — > civ) could 
itself be measured with the required precision, then it could be used to construct R2 from 
the experimentally more accessible ratio R3. 



III. RADIATIVE CORRECTIONS AT NEXT-TO-LEADING ORDER 



In Ref. we included two subsets of radiative corrections. First, we used an effective 
Lagrangian evaluated at the scale \i = rrib, 

AG F 



V2 



V, 



it b 



CMu a ^^P L b a s' l3 ^P L cp + C 2 {lj)u a ^^P L b l3 s'^P L c a +h.c. . (3.1) 
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Employing such a Lagrangian has the effect of summing all logarithms of the form 
a™ \ia n (m w / m b ) . At leading log order, Ci(m b ) ~ 1.11 and C 2 (m;,) — —0.24 f5[. Second, 
we performed a partial one-loop calculation of the radiative correction to the decay rate it- 
self. Although we included the largest contributions, our calculation was incomplete because 
we omitted terms which were not proportional to (C\ + C2/N C ) 2 . This allowed us to consider 
only gluon exchanges within color-singlet currents, simplifying the result enormously. At 
leading log order terms of this form contribute 96% of the total decay rate, so we hoped that 
the error from this approximation would not be too large. The combined radiative correction 
which we found, keeping all terms of order a™ \n n {my//m b ) and of order a s but only some 
of order a 2 \n(m\y / rrib) , was x[l + g(x c )a s /Tc], where \ came from QCD running between 
My/ and m b and g(x c ) was a finite radiative correction to the decay rate. For x c = 0.09 and 
a s = a s (mb) = 0.21, this gave 1.32, or a combined correction of about 30%. 

The largest ambiguity in this result comes from the scale fi at which the one-loop cor- 
rection to the decay process is to be evaluated. This ambiguity can be removed only by 
performing a full calculation at next-to-leading log order, including consistently all terms of 
order In" (mw/mb)- We present the results of such a calculation here. We have followed 
closely the analogous calculation of b — > cud' by Bagan et al. , and have used their partial 
results where appropriate. 

We refer the reader to the excellent exposition of Ref. || for a detailed discussion of the 
method of the analysis, and here present only our results. A brief synopsis of our calculation 
is found in the Appendix. We write the answer in the form 

Ri = \V ub /V cb \ 2 {l + n(/z, x c ) + . . .} , 

52 = \V ub /V cb \ 2 {l+r 2 (ji,x e ) + ...} , (3.2) 

where rj(/x, x c ) are of order a s . Here \x is the renormalization scale and x c = m 2 c /m 2 is the 
ratio of the heavy quark pole masses. We take m q = for q = u,d, s, which is in this process 
an excellent approximation even for the strange quark. 

Taking the reference values m b = 4.80 GeV, m c = 1.45 GeV (so x c = 0.09), and \i = 
4.8 GeV, we find 

n(fi,x e ) = 0.40 and r 2 (/x, x c ) = 0.12 . (3.3) 

Because the leading logarithms cancel in the ratio 62, r 2 starts at order a s rather than at 
order a s \a(My//m b ) and is considerably smaller than r%. In Fig. [I], we display the variation 
of ti with x c and \i. In Fig. [IJa, we vary m b between 4.5 GeV and 5.1 GeV, fixing m c by 
the heavy quark symmetry constraint m& — m c = 3.35 GeV and taking // = m b . We see 
that the dependence on x c is very mild; over the conservative range considered, r\ varies 
by ±0.01 and r 2 by ±0.02. In Fig. |l]b, we fix m b = 4.8 GeV and vary fi between \m b and 
Im.6. For r 2 , the dependence is soft, approximately ±0.02. However, for r% a significant in- 
dependence remains even at next-to-leading order. For /i as low as \m b , we have r\ = 0.53. 
We will choose to assign an asymmetrical error of (±0.10, —0.05) to the /j-dependence of r\. 
Combining the variation in x c and /1, then, we find the results 

n = 0.40±g:Jg and r 2 = 0.12 ± 0.03 . (3.4) 
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(a) 



4.0 5.0 6.0 

/j (GeV) 
(b) 



FIG. 1. Variation of with x c and [i. (a) ri(x c ) for 4.5 GeV < m& < 5.1 GeV and \i = nib. 
(b) ^(/z) for 2.4 GeV < n, < 7.2 GeV and m b = 4.8 GeV 



In the partial calculation of our previous paper, we found x(l + ga s /7r) = 1.32, to be 
compared with 1 + r\ here. We now see that this approximation underestimated the cor- 
rect next-to- leading order result by 0.08, or 20%. While our incomplete treatment gave a 
reasonable result, including the full calculation at this order turns out to be important. 



IV. BLM CORRECTIONS 

At the next order in a s , a consistent leading-log calculation requires the three- loop 
anomalous dimensions of the operators in £ and the two-loop matrix elements. However, 
since the effects of the running are not large, a useful estimate of these corrections is obtained 
by simply taking the two-loop matrix element of the singlet operator; this corresponds to 
neglecting terms of order ln(miy/m&) relative to a 2 s . A further simplification is obtained 
by only retaining the so-called "BLM" two-loop corrections [0], which are enhanced by a 
factor of /3q = 11 — 2Nf/3, where Nf = 3 is the number of light quark flavors.[| This class 
of two-loop corrections is computed easily by performing a weighted integral over the one- 
loop result calculated with a gluon mass [[!]]. While the BLM corrections are not formally 
dominant in any limit of QCD, in many processes they are found empirically to be the 
largest part of the two-loop term. In this section we calculate the BLM corrections to Ri 
and 62- 



Note that since the size of the BLM correction has nothing to do with the scale \x we used in 
the previous section to evaluate the coefficients in the effective Lagrangian, we only interpret the 
BLM correction as an estimate of the full two- loop matrix elements, not as providing information 
on the scale at which the one-loop corrections should be evaluated. We also do not include charm 
among the light quarks at the scale [i = m^. 
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These ratios require the BLM corrections to b — > ucs', b — > cud' and b — > c/z/. The 
calculation of the BLM corrections to 6 — > wcs' is identical to that for b — » ccs' || with 
one of the charm quark masses set to zero, so we refer the reader to Ref. || for details. 
The corrections are particularly simple because the corrections to the bu and cs' currents 
factorize. This feature also allows the BLM corrections to b — > cud' to be extracted easily 
from existing calculations of the semileptonic decay rate and R(e + e~ — > hadrons). Finally, 
the BLM corrections to semileptonic b —> c decays were calculated in Ref. M. Writing each 
decay rate as 

tvj, v\ i i x a sijn h ) x (a s (m b )\ 2 

T(b — > X ) oc 1 + a\ h p + . . . , (4.1) 

7T V 7T J 

we plot the af's in Fig. || for each of the relevant decays. 




0.07 0.08 0.09 0.10 0.11 0.12 0.07 0.08 0.09 0.10 0.11 0.12 



X c *c 

(a) (b) 

FIG. 2. (a) The one-loop coefficient and (b) the BLM-enhanced two-loop coefficient to 
(i) b — * ucs, (ii) b — * cud and (in) b — > civ decays. 



For the purpose of evaluating the quality of the perturbation series for the matrix ele- 
ments, one should compare this two loop result to the one loop correction g(x c ) defined in 
Eq. (|2.3| ). The reason is that neither the BLM correction nor g(x c ) has a logarithmic depen- 
dence on Mw- Hence we neglect for the moment the full NLO corrections of the previous 
section and write 



Ri = x 




+ fi(x c )Po 



7T 



7T 



+ . . . 



ft so ( a s{m b )\ 



(4.2) 



The coefficients gi(x c ) and fi( plotted in Fig. [| Taking x c 

find for Ri the results 



0.09 and \i = m b , we 
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0.07 0.08 0.09 0.10 0.11 0.12 



FIG. 3. The one- loop coefficients gi (solid) and BLM-enhanced two-loop coefficients /j (dashed). 
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71 



and for <5 2 , 



32(0.09) 



a s (m 6 ) 



2.0 



7T 



/ 2 (0.09)A) 



'a s ( m 6) 



7T 



2.1 A) 



a s (m fe ) 



7T 



1.4 A) 



0.21, 

f g s (m 6 ) 

0.14, 

f g 5 (m & ) 

I 7T 



0.09 



(4.3) 



0.06 



(4.4) 



While it is not formally consistent to include these corrections in the NLO calculations of 
the previous section, we can use them to shift the central values of the r^'s, 

n(/i = m b , x c = 0.09) -»• n + 1.09 x /i(0.09) = 0.40 + 0.10 = 0.50 , 

r 2 (fi = m b , x c = 0.09) -> r 2 + / 2 (0.09) = 0.12 + 0.06 = 0.18 . (4.5) 

Note that both corrections are somewhat larger than the error estimates from varying the 
renormalization scale in the previous section. Varying m b between 4.5 GeV and 5.1 GeV 
yields an additional uncertainty, which we estimate to be (+0.005, —0.025) on t\ and ±0.015 
on r 2 . 



V. NONPERTURBATIVE CORRECTIONS 

In addition to the perturbative corrections discussed so far, there are also nonperturbative 
contributions to R\ and 5 2 which are sensitive to the configuration of the initial B meson. 



As discussed in Ref . |L] , the cancellation of the tree level phase space factor eliminates terms 
in Ri of order Aqci)/ m b- Furthermore, while there could be in principle subleading terms at 
order l/m| proportional to the HQET parameters Ai and A2 ||10|| , these cancel in the ratio 
as well. 

This can be understood by examining the possible sources of dependence on the charm 
quark mass. These corrections enter the calculations of b — > ucs' and b — > cud' either from 
kinematic phase space functions or from performing spin sums. Consider the Fierzed form 
of the effective Lagrangian 



AG 



C -~^ Vuh 



Ci(p)sfp'Y ll P L b a UcTfPLCp + C 2 (/i)s^7 M P L 6 /3 s'^u^P^ + h.c. . (5.1) 



Substituting s' — > d! and V u b — > V c b, and exchanging u and c, we obtain the Lagrangian 
responsible for b — > cud'. Recall that we take = m s = 0. The calculation of the decay 
b — > cud' involves the computation of the polarization tensor with a massive c quark and 
a massless u antiquark, whereas the b —>■ ucs' decay involves a massless u and a massive 
c. Hence the Fierz transformation by itself does not guarantee the cancellation of all terms 
in the ratios Ri. However, because of the VV — A A structure of the polarization tensor, 
and the fact that VV and AA correlation functions are symmetric under m u <-> m c ||11|| , the 
total decay rate is also symmetric under interchange of the quark masses. Thus, both the 
phase space and l/mf corrections cancel in <5 2 . To extend the argument to Ri, note that in 
the Fierzed form and at tree level, the decay b — > cud' is the same as b — > utv T) for which 



the nonperturbative corrections were computed in Ref. |L2]. This discussion elaborates the 



observation made in Ref. n that these corrections cancel. A general argument that the terms 



proportional to Ai (but not to A2) cancel in all ratios of B decays was given in Ref. [13 . 

There remain, however, mixed terms proportional to a s A 2 /7«b, which need not cancel 
in the ratios. Recall that A2 = (B\b v ga ■ Gb v \B) /12Mb, where b v is the effective b field 
of HQET, is related to the hyperfine interaction of the b quark chromomagnetic moment 
with the light degrees of freedom in the B meson |I(J[n|. The chromomagnetic operator is 
obtained by attaching a gluon to one of the quarks in the final state, before the operator 
product expansion is performed. In the effective theory defined by the Lagrangian fl3.1|) , 
terms of order a s A 2 come from two sources. First, they arise from one-loop radiative correc- 
tions to the operator product expansion itself, in which case they are quite small. Second, 
the color structure allows terms proportional to C1C2 to arise at tree level from the inter- 
ference of the color singlet and color exchanged operators. These terms are really of order 
a s \n{My/ / 771^X2/ ml and hence are enhanced over the others. They were first calculated for 
the decay b — > cud' in Ref. \T3\. 



In fact, a straightforward argument shows that these terms are equal in size and of the 
opposite sign in b — > ucs' as compared to b — > cud' . Consider again the Fierzed form of the 
effective Lagrangian Because of the color structure when the operators are written 

in this form, the term proportional to CiC 2 A 2 is generated by attaching a gluon to the 
cu loop, performing the operator product expansion, and extracting the chromomagnetic 
moment operator ba^G^b. This term is odd under the exchange u — > c and c — >• u, 
or equivalently m\ «-> m 2 cl which can be seen immediately by inspection of the relevant 
Feynman diagrams. Alternatively, simply note that the quark and antiquark produced by 
the left-handed current carry opposite magnetic moments. To obtain the result for b — > ucs', 
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we take the limit m u = 0; as noted above, for b — > cud', we can take m c = and m u — > m c . 
Since the general result is odd in (m 2 — ml), the two limiting results are the negatives of 
each other, as promised. 

We write the result as fractional corrections to R\ and 8 2 (suppressing for the moment 
the radiative corrections of the previous sections), 

Ri = \V ub /V cb \ 2 {l + £ 2 (fi,x c ) + ...}, 

S 2 = \V ub /V cb \ 2 {l + 2£ 2 (fx,x c ) + ...}, (5.2) 

where 

e 2 (n, Xc ) = cmcm m ~* c)3 ^ (5.3) 

and f(x c ) = 1 — 8x c + 8x% — x\ — 12x 2 \yix c is the tree level phase space function. For 
consistency, C\{p) and C 2 {\i) should be evaluated at leading order, to include only terms of 
order a™ ln n (M w /m b )\ 2 /m 2 . The scale dependence of A 2 is given by JL4[] 

A 2 (/i) = \ 2 {m h ) \a s ([i)/a s (m b )] 1 . (5.4) 



For /i = m b = 4.8 GeV, and with X 2 (m b ) = 0.12 GeV fixed by the B — B* mass splitting |T4 



we find £ 2 = —0.036. The variation of £ 2 with x c and /i is shown in Fig. f|. We see that, 



-0.03 



-0.04 



-0.05 
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(a) 



4.0 5.0 6.0 

H (GeV) 
(b) 



FIG. 4. Variation of t 2 with x c and [i. (a) £ 2 (x c ) for 4.5 GeV < m b < 5.1 GeV and fi = m b . 
(b) £ 2 (n) for 2.4 GeV < n < 7.2 GeV and m b = 4.8 GeV. 

as with the radiative corrections, the only significant uncertainty comes from the choice of 
renormalization scale. We assign the tiny error ±0.001 to £ 2 from the variation with x c , and 
the larger asymmetrical error (+0.010, —0.015) to the variation with fi. Combining these in 
quadrature, we see that the first is negligible, and we estimate 



-0.036; 



-0.010 

-0.015 • 



(5.5) 
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Noting that the fractional correction to 82 is 2£ 2 and comparing to the radiative correction 
r 2 (|3.4| ), we see that this is actually an important uncertainty for <5 2 . 

The leading purely nonperturbative corrections to R x and <5 2 arise at order l/mf. The 
largest such corrections are associated with annihilation processes such as bu — * cs', since 
they are enhanced by a relative phase space factor of 16ir 2 . They are also spectator depen- 
dent, contributing differently to B~ and B° decays. In Ref. fll], we discussed the derivation 
and computed the annihilation terms in R x . Here we will recall those results, as well as 
present results for <5 2 . As before, we will present a fractional correction, of the form 

Ri = \V ub /V cb \ 2 {l + ai (B-,B°) + ...}, 

82 = \V ub /V cb \ 2 {l + a 2 (B-,B°) + ...}, (5.6) 

where all other corrections have been momentarily suppressed. The terms at depend on 
nonperturbative matrix elements of four-quark operators, parameterized by "bag factors" 
B % and ti JT§: 

(B q \ h'fqL qLlpbh \B q ) = \ P B m 2 Bq B x , 
(B q \ b R q L q L b R \B q ) = \ f 2 B m 2 Bq B 2 , 
(B q \ b L ^T a q L q Ll ,T a b L \B q ) = \ f 2 m 2 Bq e x , 

(B q \ b R T a q L q L T a b R \B q ) = \ f 2 m 2 Bq e 2 . (5.7) 

In the vacuum insertion ansatz, only color single operators contribute to the decay and we 
have B\ = B2 = 1 and ei = e 2 = 0. More generally, the color octet parameters q are of 
order 1/N C in the limit iV c — > 00. 

In terms of these parameters, we find the corrections 

ai(fl-) = 16?r2/ f ^ " ^ {(C x + \C 2 ) 2 [{1 + 2x c )B 2 - (1 + \x c )B x ] 

m b J\ x c) 

+ |C 2 2 [(l + 2a; c )e 2 -(l + ix c )e 1 ]}, 
a x (B°) = 16?r2/ f S ~ x ^ sin 2 6c Uc 2 + |d) 2 [(l + 2x c )B 2 - (1 + \x c )B x \ 

m b J\ x c) 

+ §C 2 [(1 + 2a; c )e 2 - (1 + \x c )e x ]} (5.8) 

for R x , and 

« 2 (S-) = J f;, cJ {(d + |C 2 ) 2 [(1 + 2x c )5 2 - (1 + \x c )B x ] 

m b J\ x c) 

+ |C 2 [(l + 2x c )e 2 -(l + ix c ) ei ] 

-|(C? + 6C X C 2 + C 2 )B X - 2(C 2 + C|)ei} , 

1 fi?r 2 i 2 (\ — r \ 2 

a 2 (B°) = - Jb 2 ) Xc) cos 2^ {(C 2 + \C x f[{l + 2z c )£ 2 - (1 + l^)^] 

m b J\ X c) 

+ |C?[(1 + 2x c )e 2 - (1 + \x c )e x )} (5.9) 
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for 82- In deriving the corrections to the denominator of <5 2 , we have adapted the results 
of Ref. H for the channel b -> cud'. With m b = // = 4.8 GeV, / B = 200 MeV and 



sin9c = 0.22, we find 

oi(S~) = -0.48Si + 0.54^2 - 0.018ei + 0.021e 2 
Oi(B°) = -0.00034^! + 0.000385 2 - 0.018ei + 0.020e 2 
a 2 {B-) = -0.43.Bi + 0.54^2 - 1.146! + 0.021e 2 , 

a 2 (B°) = 0.00635! - 0.0071E 2 + 0.34ei - 0.38e 2 . (5.10) 

Note that in a 2 (B^, B°), the coefficients of e« are quite large. This reflects the potentially 
significant contribution of color-octet annihilation processes to nonleptonic B decays, as 
observed in Ref. ||15|| . Since the parameters 6j are not known well, the large size of these 
terms introduces a problematic uncertainty into the denominator of <5 2 . If we take the 
vacuum insertion ansatz, in which 6j do not contribute, we have 

a x {B-) = 0.062, at(B°) = 4.4 x 10" 5 , 

a 2 {B-) = 0.112 , a 2 (B°) = -8.1 x 10~ 4 . (5.11) 

Unfortunately, it is hard to assess the uncertainty due to nonzero e,. For want of a better 
procedure, let us survey briefly the available models for estimating the matrix elements ( p.7| ). 
The most reliable of these, in principle, is the lattice QCD result JH 



£i(m 6 ) = 1.06 ± 0.08 , ei(m 6 ) = -0.01 ± 0.03 , 

B 2 {m b ) = 1.01 ± 0.06 , e 2 (m ft ) = -0.02 ±0.02, (5.12) 

where the quoted errors include neither quenching errors nor the systematic uncertainty due 
to the extrapolation to the chiral limit. Both of these issues can be addressed in future, 
more precise calculations. There also exist calculations in the framework of QCD sum rules, 



which give 17 



and 18 



fii~l, ei ~ -0.15, 

B 2 ~l, e 2 ~0, (5.13) 



5i(m 6 ) = 1.01 ± 0.01 , ea(m 6 ) = -0.08 ± 0.02 , 

B 2 {m b ) = 0.99 ± 0.01 , e 2 (m fe ) = -0.01 ± 0.03 , (5.14) 



as well as an HQET QCD sum rule calculation which yields [19 

5i(m 6 ) = 0.96 ± 0.04 , ei(m 6 ) = -0.14 ± 0.01 , 

B 2 {m b ) = 0.95 ±0.02, e 2 (m fe ) = -0.08 ± 0.01 . (5.15) 

While the 5j's are consistently within 5% or so of unity, there is a large spread in the values 
of the 6j parameters. Of course, it makes no sense to "average over models" as a method for 
assigning values to Bi and 6j. Instead, we adopt the procedure of taking the lattice QCD 
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results as central values but inflating the errors both to be conservative and to reflect the 
variety of values which QCD sum rules yield for 6j. To be even more conservative, we inflate 
the errors symmetrically, so that we use the sum rules results to set the magnitude, but not 
the sign, of the uncertainty in the lattice calculations. The central values and errors which 
we choose are then 

#i(m 6 ) = 1.06 ± 0.10 , ei(m 6 ) = -0.01 ± 0.10 , 

B 2 (m b ) = 1.01 ±0.10, e 2 (m b ) = -0.02 ±0.10, (5.16) 

One could imagine a less conservative procedure, especially if one had particular confidence 
in one of the calculations quoted above, but this is not the approach which we will follow. 

Inserting these parameters ( |5.16| ) into the solution ( |5.10| ) for the nonperturbative correc- 
tions, we find 

a x (B-) = 0.04 ± 0.07 , a^W) = ± 0.003 , 

a 2 (B-) = 0.10 ±0.13, a 2 (B°) = ± 0.05 . (5.17) 

Note that the fractional error associated with the weak annihilation contribution is partic- 
ularly large for S 2 in charged B decays, due to its enhanced sensitivity to e^. 



VI. PHENOMENOLOGY OF V ub FROM i?i AND R 2 



Combining the results of the previous sections, we now present estimates for central 
values and uncertainties for Ri and 5 2 = 1 — R 2 . Due to the sizable flavor-dependent 
corrections associated with the spectator contributions, we present our results separately for 
charged and neutral B decays, which we denote by introducing a suitable superscript. Our 
results take the form 



s 2 

5°, 



\V ub /V cb \ 2 1 + r 1 (x c ,fi) +£ 2 (x c ,ii) + a x (B ) 
\V ub /V cb \ 2 \l + rx(x c , fx) + £ 2 {x c , /i) + a,i(B ) 
\V ub /V cb \ 2 1 + r 2 (x c , n) + 2£ 2 (x c , fi) + a 2 (B~ 
\V ub /V cb \ 2 \l + r 2 (x c , n) + 2£ 2 (x c , //) + a 2 (B° 



(6.1) 



where comes from perturbative QCD radiative corrections and £ 2 and a, represent frac- 
tional corrections due to nonperturbative effects. Putting our results together, we find 





\v ub /v cb 


2 


"i cn+ - 10 + - 005 +°- 010 

i - JU -0. 05-0.025-0. 015 


*? = 


\v ub /v cb 


2 


1 /|«+0. 10+0.005+0.010 
i -^ :O -0. 05-0.025-0.015 


5 2 = 


\V ub /V cb 


2 


1.21 ±0.03 ±0.015 


5° - 


\v ub /v cb 


2 


1.11 ±0.03 ±0.015 



±0.005 



+0.015 + n i o 
-0.030 ^ U - ±D 



+0.015 
-0.030 



±0.05 



(6.2) 



In these expressions, the first error is our estimate of the uncertainty from NLO pertur- 
bative QCD corrections, second is due to uncertainty in the BLM part of the two-loop 
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corrections, the third represents uncertainty in the a s \\i{Mw jm^^lm^ term, and fourth 
is due to spectator-dependent 1/mf effects. We expect the net effect of other l/mf effects, 
not enhanced by the phase space factor of 16n 2 , to be safely below the level of the estimated 
uncertainty. 

We can combine the errors quoted in Eq. ( |6.2| ) by taking into account the correlations 
among the various sources of theoretical uncertainty. However, we would like to emphasize 
that this procedure of estimating and combining theoretical errors, while widespread, is 
purely conventional and has no rigorous statistical meaning. With this in mind, our best 
estimate of the central values and overall uncertainties is 



\v ub /v c 



ct>\ 



\V u b/V c 



cb\ 



1.50 ±0.15 
1.46 ±0.10 



6; 



\v ub /v c 



cb\ 



5° 2 = \V ub /V cb \ 



1.21 ±0.15 
1.11 ±0.10 



(6.3) 



Although the uncertainties are generally smaller for the neutral B decays, in all cases 
the theoretical errors are at approximately the level of ten percent. For Ri, the error 
is dominated by residual uncertainties in the next-to-leading order radiative corrections. 
Reducing them substantially would require no less than a next-to-next-to-leading order 
calculation. For <5 2 , the uncertainties come primarily from the poorly known strong matrix 
elements needed for the annihilation contributions, especially from the color octet bag factors 
ej. The best prospect for improvement here is in a future generation of unquenched lattice 
calculations. Were these to become available, overall theoretical errors in 5 2 at the five 
percent level would be within reach. 

In summary, we have studied the possibility of extracting the CKM matrix element V ub 
from ratios of inclusive nonleptonic B decays. We have shown that the ratios of inclusive 
decay widths defined by Eqs. ( |2.1|) and (|2.4j ) are impressively "clean" theoretically. We 
have estimated the impact and uncertainty associated with the NLO radiative corrections, 
and have included the BLM part of the two-loop term. In addition, we have studied the 
impact of the leading non-perturbative corrections on Ri and R 2 . There is no doubt that 
the measurement of the Ri and R2 would be a challenging enterprise. Unfortunately, the 
somewhat experimentally easier ratio of R3 of Eq. (|2.6|) has larger theoretical uncertainties, 
from radiative corrections which would need to be computed at better than the 1% level 
before the method could be used to extract V ub . Nonetheless, the ratios Ri and R 2 of 
nonleptonic decay widths offer a new and tantalizing approach to measuring the important 
but poorly known CKM matrix element V ub . 
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APPENDIX A: 

Here we outline the procedure for calculating radiative corrections to Ri and 82 with next- 
to-leading log (NLO) accuracy. This calculation amounts to the computation of perturbative 
corrections to b — > ucs 1 . We follow closely the procedure outlined in Ref. |J. 
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FIG. 5. Feynman diagrams for the calculation of QCD corrections to b 
represent insertions of the operator 0\ or O2. Dashed lines represent gluons. 



ucs' . Black dots 



A full NLO calculation of QCD radiative corrections to b — > ucs' involves the compu- 
tation of the eleven loop diagrams depicted in Fig. [5|. These diagrams can be conveniently 
"packaged" into five classes. In fact, only the diagrams of a single class will have to be 
computed, as the others may be extracted from existing calculations of perturbative QCD 
corrections to polarization tensors and inclusive semileptonic b decays. As discussed in 
Ref. 0, this is achieved by the application of Fierz relations to 0\ = Uori^P^bas'^P^c^ 
and O2 = Ua^fiPibp s'^PlCo. (note that our choice of 0\ and O2 is opposite to that of 
Ref. ||). Although in general Fierz symmetry is broken by regularization, Naive Dimen- 
sional Regularization (NDR) along with the choice of evanescent operators advocated in 
Ref. E0| actually preserves Fierz relations for renormalized operators. This allows us to use 
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the Fierz transformation on a diagram by diagram basis 0, which significantly simplifies 
the task of computing the graphs of Fig. [| in the limit of massless u and s quarks. 

We use the effective Lagrangian of Eq. ( p.lQ calculated at NLO accuracy. It is convenient 
to define the Wilson coefficients C±(fi) of the multiplicatively renormalized operators 0± = 
\{Oi ± O2). At this order, they are given by 



C±(Ai) = Ci0i)±C 2 (Ai)=L±0i) 



1 + 



4vr 2/5 
a s (m w ) - a s (jj) 



a s {m w ) - a s (fj) 7± I 7± A 

' T f /9b 



H ; -d± 



47T 



4tt 



47T 



(Al) 



where for Nf flavors and N c = 3 colors the anomalous dimensions 7 ± = 7 ± (a s /47r) + 
7i 1) (a s /47r) 2 + ... are ^ ] = 4, 7 i 0) = -8, 7^ = -7 + 4AT//9, and 7 i 1} = -U-8N f /9. The 
QCD (3 function is given by 



P 



+ 



(3 = n--N f , 



h = 102 - y JV) . 



(A2) 



The coefficients of Oj- at leading logarithmic order are 

7i 0) /2A, 



Os(%)\ 



(A3) 



Finally, the result of two loop matching at fi = for the effective Lagrangian (in NDR) 
is contained in 



±B 



2N r 



where 



As advocated in Ref. ||, it is useful to combine Eq. 

R± = B± + 



(0) / (l) a 
7± / 7± Pi 



2Po \ 7 i 0) Po 



B = 11. (A4) 
with the anomalous dimensions, 

(A5) 



so that R± are independent of the renormalization scheme. Following Ref. ||, the decay 
rate for b — > ucs' can be expressed as 



r(6 -> ucs') 



G 2 F ml\V ub \ 2 

1927T 3 



2L 2 ,( / i) + L 2 » 



TV 



2L 2 h (/i) J R+ + L 2 (//)#_ 



+ 
+ 



2tt 
3vr 



[L+(/i) + ^_(/i)] 2 c 22 (a; c ) + [£+(//) - L_(//)] 2 c u (x c ) 
L+(/i) 2 - L 2 (ji)] c 12 (x c ) \ , 



(A6) 
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G a( x c) 
G e (x c ) 



0.07 0.08 0.09 0.10 0.11 0.12 

X c 

FIG. 6. Variation of Gs a ^ e \ with x c . The dashed line is G a (x c ), the dash-dotted line is Gb(x c ), 
and the solid line is G e (x c ). 



where we define cn(x c ) = G c + Gd, C22(x c ) = G a + Gb, and ci 2 (x c 
Here the five classes of graphs Gi are defined as 



G a + Gb + G e 



B. 



7T 
7T 

7T 



7T 

7T 



K Im 
K Im 
K Im 
K Im 
K Im 



(2) + (2)t + (3) + (9) + (9)t + (12) 
(4) + (5) + (7); , 

(2) + (2)t + (5) + (ll) + (ll)t + (12) 

(3) + (4) + (6); , 

(4) + (8) + (10) + (10) t + (ll) + (ll) t 



(A7) 



with-K" = 1927r 3 /m^/(x c ) and f(x c ) = 1— 8x c +8xl~ x\ — 12x 2 In x c . Examination of Eq. (|A7|) 
reveals that G a = G c and can be extracted from the calculation of QCD corrections to the 

0, 



semileptonic decay b 

G a {x c ) = — {4(1 - x c )(75 - 539x c - 476a; 2 + 18a; 



48/(x c 



- 16tt 2 (3 - 24x c - 36a; 2 + 16x 3 c - 2x A c ) - 3456a; 2 (((3) - Li 3 (z t 

- 96(1 - 8x c + 36a; 2 + 16x 3 c - 2x A c ) U 2 {x c ) 

- 8(1 - x 2 )(31 - 320x c + 31x 2 ) ln(l - x c ) 

- A8(2x c + 15a; 2 - f x 3 c + fx 4 c - 8n 2 x 2 c + 24x 2 Li 2 ( 



[Xr 



+ 2(1 -x 2 )(l-8x c + x 2 )ln(l 



lnx c | . 



(A8) 



Likewise, Gb = Gd and can be extracted from the computation of perturbative QCD 
corrections to the polarization operator of vector and axial currents [p| , p^| , 



°b( x c) = AQ l J (l " x c )(18 - 476a; c - 539a; 2 + 75^ 
48/(x c ) L v 
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+ 4x 2 (36 + 8x c - x 2 c )(n 2 - 3 In 2 x c ) 

- 2(1 - x\) (31 - 320x c + 31x 2 c - 12(1 - 8x c + x\) lnx c ) ln(l - x c ) 

- 2x c (132 + 90a; c - 308x 2 + 31x^) lnx c 

+ 24(2 - \§x c - 36x 2 + - x\ + 12a; 2 lnx c ) Li 2 (x c ) 

+ 864x 2 (c(3)-Li 3 (x c ))}. (A9) 

Hence G e is the only new quantity which needs to be computed. We have done this numer- 
ically, modifying a subroutine given to us by P. Ball ||. Since each of G a ,Gb,G c , and Gd 
is related to a physical process, the individual classes of the diagrams are gauge invariant 
and free of infrared singularities. We present plots of Gi(x c ) in Fig. ||, with /i = m h and 
m b — m c = 3.35 GeV. We see that the dependence on the quark masses is not very strong. 
To obtain the radiative correction to Rx, the results of this calculation must be combined 



with the one- loop result for b — > civ [23]. For the denominator of 62 , we use the result of 



Ref. 0. The radiative corrections to the two observables are plotted in Fig. |T] of the text. 
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